Analytical model of the enhanced light transmission through subwavelength metal 
slits: Green's function formalism versus Rayleigh's expansion 
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We present an analytical model of the resonantly enhanced transmission of light through a sub- 
wavelength nm-size slit in a thick metal film. The simple formulae for the transmitted electromag- 
netic fields and the transmission coefficient are derived by using the narrow-slit approximation and 
the Green's function formalism for the solution of Maxwell's equations. The resonance wavelengths 
are in agreement with the semi-analytical model [Y. Takakura, Phys. Rev. Lett. 86, 5601 (2001)], 
which solves the wave equations by using the Rayleigh field expansion. Our formulae, however, show 
great resonant enhancement of a transmitted wave, while the Rayleigh expansion model predicts 
attenuation. The difference is attributed to the near-field subwavelength diffraction, which is not 
considered by the Rayleigh-like expansion models. 

PACS numbers: 42.25.Bs, 42.65.Fx, 42.79.Ag, 42.79.Dj 



I. INTRODUCTION 

The near-field localization and resonantly enhanced 
transmission of light by metal subwavelength nanos- 
tructures, such as a single aperture, a grating of aper- 
tures and an aperture surrounded by grooves, attract 
increasing interest of researchers. The recent studies 
QIIEIEEIIILIEIEIEIE EH have pointed out that the 
localization and enhanced (anomalous) transmission of 
light by a grating of holes or slits can be better un- 
derstood by elucidating the optical properties of a sin- 
gle subwavelength slit. Along this direction, it was al- 
ready demonstrated that a TM-polarized light wave can 
be localized in the near-field subwavelength zone of a 
metal slit and simultaneously enhanced by a factor of 

about io-io 3 HlSEESHSIISmHm. 

The general analysis and interpretation of these results, 
however, are very complicated because the studies are 
based mainly on purely numerical computer models. Re- 
cently, a simple semi-analytical model of the light trans- 
mission by a subwavelength metal slit was developed 
0. The study clearly showed that the transmission co- 
efficient versus wavelength possesses a Fabry-Perot-likc 
resonant behavior. Unfortunately, the model [l| pre- 
dicts transmission peaks with very low magnitude (at- 
tenuation) jwhile the experiments and computer models 
0, El EL IE IE IE IE E3 demonstrate the great resonant 
enhancement of a transmitted wave. 

In this paper, we present an analytical model of the 
resonantly enhanced transmission of light through a sub- 
wavelength nm-size metal slit. The simple formulae for 
the transmitted electromagnetic fields and the transmis- 
sion coefficient arc derived by using the narrow-slit ap- 
proximation and the Green's function formalism for the 
solution of Maxwell's equations. The article is organized 



as follows. The model and formulae are described in sec- 
tion [n] In section ITTT1 the model is compared with the 
results of the semi-analytical model 0. We show that 
the resonance wavelengths determined by our formulae 
are in agreement with the model pj, which solves the 
wave equations by using the Rayleigh field expansion. 
The formulae, however, indicate great resonant enhance- 
ment of a transmitted wave, while the Rayleigh expan- 
sion model predicts attenuation. The summary and 
conclusion are given in section llVl 



II. ANALYTICAL MODEL BASED ON 
GREEN'S FUNCTION FORMALISM 

a. Theoretical background Let us consider the scat- 
tering of TM-polarized light by a subwavelength slit in a 
thick metallic film of perfect conductivity. From the lat- 
ter metal property it follows that surface plasmons do not 
exist in the film. Such a metal is described by the classic 
Drude model for which the plasmon frequency tends to- 
wards infinity. We follow the Neerhoff and Mur theoreti- 
cal development based on the Green's function formalism 
for the solution of Maxwell's equations [111 Il2l 1 1 3| . 

The schematic diagram of the light scattering is shown 
in Fig. ^ The transmission of a plane wave through a 
subwavelength (A > 2a) slit of width 2a in a metal film of 
thickness b is considered. In region I, the incident wave 
propagates in the (x—z) plane at an angle O with respect 
to the z axis. The magnetic field of the wave is assumed 
to be time harmonic with the frequency to — kc and both 
polarized and constant in the y direction: 

H(x, y, z, t) = U(x, z)exp(—iujt)e y . (1) 
The electric field of the wave is found from the scalar field 
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FIG. 1: Schematic diagram of light scattering by a subwave- 
length slit (waveguide) in a thick metallic film. 



U(x, z) using Maxwell's equations. Thus, the electromag- 
netic field is determined by a single scalar field U(x,z). 
In the regions I, II and III, the field is represented by 
Ui(x, z), U 2 {x,z) and U 3 (x, z), respectively. The field 
satisfies the Helmholtz equation (V 2 + k 2 )Uj = 0, where 
j = 1,2,3. In region I, the field U\(x,z) is decomposed 
into three components, U\(x, z) — U l (x,z) + U r (x,z) + 
U d (x, z), each of which satisfies the Helmholtz equation. 
The incident field U l (x,z) — exp[i/ci(a;sin0 — zcosG)] is 
assumed to be a plane wave of unit amplitude. U r (x, z) = 
U l {x, 2b — z) denotes the field that would be reflected if 
there were no slit in the film. U d describes the field 
diffracted by the slit into region I. To find the field, the 
2-dimcnsional Green's theorem is applied with one func- 
tion given by U(x,z) and the other by a conventional 
Green's function: 



(V 2 + fc 2 )Gj = -S(x -x',z- z'), 



(2) 



where (x, z) refers to a field point of interest; x' and z' 
are integration variables, j = 1,2,3. Since Uj satisfies 
the Helmholtz equation, Green's theorem reduces to 



U d {x,z) 
U 3 (x, z) 
U 2 (x, z) 



-— / Gi(x,z;x',b)DU b (x')dx', (4) 

£2 J-a 

— [ G 3 (x,z;x',0)DU (x')dx', (5) 

e 2 J-a 



[G 2 (x,z;x',0)DU (x') 
-U (x')d z >G 2 (x, z; x\ z')\ z - t0 +]dx' 
+ f [G 2 {x,zyx',b)DU b {x') 

J —a 

-U b (x')d z >G 2 (x,z;x', z')\ z ^ b -]dx' , (6) 



where the boundary fields are defined by 

U Q {x) = U 2 (x,z)\ z ^ Q +, 

DUq(x) ee d z U 2 (x,z)\ z ^ +, 

U b {x) EE U 2 (x, z)\ z ^ b - , 

DU b {x) ee d z U 2 (x,z)\ z ^ b - . 



(7a) 
(7b) 
(7c) 
(7d) 



In regions I and III the Green's functions are given by 



Gi (x, z; x 1 , z') = i [H™ (hR) + (hR')] , (8) 
G 3 (x,z;x',z') = l -[H^(k 3 R)+H^(k 3 R'% (9) 



with R = [{x-x') 2 + {z-z') 2 ] l l 2 , R! = [(x-x') 2 +(z+z'- 

2b) 2 ] 1 ' 2 , and R" = [{x - x') 2 + (z + z') 2 } 1 / 2 . Here, 
is the Hankel function. Inside the waveguide (region II), 
the method of images can be used; the Green's function 
is given by the waveguide multimode (m = oo) expansion 



G 2 (x, z;x',z') 



4a7o 



exp(i7o|z - z'\) 



i 

+ 2a~ 

X cos 



7m lex P(nmN - z'\) 

m— 1 

mir(x' + a) mir{x + a) 



2a 



■ cos ■ 



2 a 



(10) 



U{x,z)= j (Gd n U - Ud n G)dS, (3) 

J Boundary 

where n is the normal vector at the boundary surface. 
The unknown fields U d (x,z), U 3 (x,z) and U 2 (x,z) re- 
spectively for the regions (b < z < oo), (— oo < z < 0) 
and (|x| < a, < z < 6) are found using the reduced 
Green's theorem and the standard boundary conditions 



where 7 m = [k 2 — (mn/2a) 2 } 1 ^ 2 . The field can be found 
once the four unknown functions in Eqs. I7al7dl have been 
determined. The functions are determined by the four 
integral equations [l3j : 

/a 
Gx(x,b; x' ,b)DU b (x')dx' , (11) 
-a 

U ( x ) = — ( G 3 (x,0;x',0)DU a (x')dx', (12) 

e 2 J-a 
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~G 2 (x,b;x',0)DU (x') 
U (x')d z iG 2 {x, b; x', z')\ z ^ +]dx' 

G 2 {x,b;x',b)DU b (x')dx', (13) 



;Uo(x) 



[G 2 (x,0;x',b)DU b (x') 
U b (x')d z ,G 2 (x,0;x' , z%^ b -]dx' 

G 2 (x,0;x',0)DU (x')dx\ (14) 



where \x\ < a, and Ul(x) — exp[ifci(a;sin8 — 6 cos 6)]. A 
set of the four coupled integral equations can be solved 
numeri cally by dividing each integral in the expressions 
(|11I14|) into jV subintervals of width 2a /N (for more de- 
tails, see |l3J). The four boundary functions are found 
through the following matrix equation: 

IIoIII TjII 



11 + Is 1 

S 11 + IS 111 D 1 ^ 1 - R 11 



D^S 111 R 

i 

2 





'DUo(x)~ 




Ul{x) 




DU b (x)_ 




_2D u Ul(x)_ 



(15) 



where 



R kj - 



2A 7( 



■ exp(«7 6) 



ly m =l 7m 
2^ exp(i 7o 6) 

^ OO 



LXJ 

N ^ lm m ' 

m— 1 



2N l0 
2N 



■ cos 



mn(j - 1/2) 



X cos 



N 

rrm(k - 1/2) 



TV 





~rmr~ 


sin 


L 2 JV" J 





(16a) 

(16b) 
(16c) 

(16d) 



b. The analytical model and formulae Let us now 
present an approximate analytical solution of the scat- 
tering problem. For a subwavelength slit, the narrow-slit 
condition 2o <C A is a very good approximation. In such 
a case, it can be easily demonstrated that an accurate 
solution exists for the above matrix equation in the case 
of N = 1, which is valid at \z\ > 2a /N (see, ref. [Ts]). We 
have derived the simple formulae for the transmitted elec- 
tromagnetic field and the transmission coefficient using 
the narrow-slit approximation. After simple calculation 
using the above mentioned conditions, the analytical for- 
mulae for the magnetic H (x, z) and electric E{x, z) fields 
in region III can be presented as 

H{x,z) = iaH^\k^/x 2 + z 2 ){DU )ie y , 



E x (x, z) = -a 
E y {x,z)=0 



\/x 2 + z- 



(17) 

-.HPiW^ + z^DUo)!, (18) 

(19) 



and 



E z (x, z) = a 



* rr(D 



\Jx 2 + 



H[ 1 \kVx 2 + z 2 ){DU )i, (20) 



where 



exp(ikb — ikb cos ( 



ik [exp(ifcb) • (A - ik)] 2 - (A + ik) 2 



(21) 



and 



A = 



H, 



'P (*a) + (|) [Fo(fca)^i 1} (ka) 

-Pi (ka) (ka)}] . (22) 



Here, H^ 1} (ka) and H[ 1} (ka) are the Hankcl functions, 
and Fo(ka) and F\(ka) are the Struve functions; for the 
sake of simplicity, we presented the formulae for the slit 
placed in vacuum. The transmission coefficient is deter- 
mined by the time averaged Poynting vector (energy flux) 
S of the electromagnetic field. The vector is calculated 
as S = (E x H* + E* x H). The transmission coefficient 
T is defined as the integrated transmitted flux S 3 (z) di- 
vided by the integrated incident flux S l . The flux S 3 (z) 
is integrated at the distance z (z > 2a) from x m i n = -co 
to x max = oo, and the flux S 7, is integrated over the slit 
width 2a at the slit entrance. For the transmission coef- 
ficient we have obtained the following simple analytical 
formula: 



T(X,a,b) 



kcosO 



[(Re(DUo)i) 2 + (Im(DUo)i) 2 ] . (23) 



Notice, that in agreement with the energy conservation 
condition, the transmission coefficient does not depend 
on the distance z. This means that the near-field and far- 
field transmission coefficients are of equal value, T(2a < 
z<\) = T(z>\). 



III. COMPARISON WITH MODEL BASED ON 
RAYLEIGH'S FIELD EXPANSION 

The analytical formulae (|17I23|) for the transmitted 
electromagnetic field and the transmission coefficient 
have been derived by using the narrow-slit (2a <C A) ap- 
proximation and the Green's function formalism for the 
solution of Maxwell's equations. The main objective of 
our analysis is to test the accuracy and range of validity 
of the formulae. We have calculated the magnetic H(x, z) 
and electric E(x, z) fields of the transmitted wave and the 
transmission coefficient T for different parameters of the 
slit. As an example, Figs. 2-5 show the numerical results 
for the incident angle equal to zero. Figure 2(a) shows 
the field distributions \H y (x,z)\, \E x (x,z)\ and \E z (x,z)\ 
calculated by the analytical formulae (|17I20I) for the slit 
width 2a — 150 nm and the screen thickness b = 150 nm. 
The field distributions calculated at the distance \z\ = 2a 
by the numerical evaluation J13J of the integral equations 
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FIG. 2: The field distributions \H y \, \E X \ and \E Z \ of the 
transmitted wave with A = 2500 nm calculated at \z\ = 2a 
for the slit width 2a — 150 nm and the screen thickness 
b = 150 nm: (a) analytical model and (b) rigorous numer- 
ical model. 



(|11I14|) are shown in Fig. 2(b) for the comparison. 
We notice that the distributions are practically undistin- 
guishable at \z\ — 2a. The calculations show that the 
accuracy of the analytical formulaes increases with fur- 
ther increasing the distance \z\. The field components 
H y and E x are collimated to approximately the aperture 
width 2a (see, Fig. 2). Thus, the energy flux S z of a 
wave passing through a subwavelength slit can be used 
to provide a subwavelength image with the optical reso- 
lution of about 2a. We calculated the magnetic H (x, z) 
and electric E(x, z) fields of the transmitted wave also 
in the region (\z\ < 2a) by using the analytical formula 
(Fig. 3(a)) and by the numerical evaluation ^| of the 
integral equations I|llll4fl (Fig- 3(b)). The results show 
that the accuracy of the formulae I|17l20fl decreases with 
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FIG. 3: The field distributions \H y \, \E X \ and \E Z \ of the 
transmitted wave with A = 2500 nm calculated at \z\ — 0.1a 
for the slit width 2a — 150 nm and the screen thickness 
b = 150 nm: (a) analytical model and (b) rigorous numer- 
ical model. 



decreasing the value \z\ in the region (|z| < 2a). 

The transmission of light by the slit is a process that 
depends on the wavelength. Due to the dispersion, the 
amplitude of a wave does change under propagation 
through the slit. This leads to attenuation or amplifi- 
cation of the wave intensity. In our analytical model, the 
transmission coefficient T(A, a, b) is described by the ana- 
lytical formula (|23|l . We now check the consistency of the 
results by comparing the transmission coefficient T(X) 
calculated by using the semi-analytical model with 
those obtained by the formula 12- il) . Analysis of Eq. 1251 
indicates the slit-wave interaction behavior that is sim- 
ilar to those of a Fabry-Perot resonator. The minimum 
thickness of the screen is required to get the waveguide 
resonance inside the slit at a given wavelength. Figure 4 
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FIG. 4: The transmission (curve A) as a function of the 
wavelength calculated by using the analytical formula I23H 
for the slit width 2a = 150 nm and the screen thickness 
b — 150 nm. The curve B shows the transmission calcu- 
lated in the study Q. The dotted vertical line corresponds 
the resonant wavelength given by the analytical formula. 



shows the transmission coefficient (curve A) as a function 
of the wavelength A calculated by using the formula |2l5)l 
for the slit width 2a = 150 nm and the screen thickness 
(b = 150 nm) smaller than the wavelength A. The trans- 
mission coefficient (curve B) calculated in the study U 
is presented in the figure for comparison. In contrast 
to the study Q, the transmission (curve A) calculated 
by our formula exhibits resonance peak, while curve B 
demonstrates no resonances. In Fig. 4, the screen thick- 
ness has been chosen in such a way as to include one 
Fabry-Perot resonant wavelength Xpp = 2b. Notice that 
the resonance peak is red shifted in comparison with the 
respective Fabry-Perot wavelength Xpp — 300 nm. In 
Fig. 5, the slit includes several Fabry-Perot resonances 
at the wavelengths Xfp — 2b/n. Figure 5 shows the 
transmission coefficient as a function of the wavelength 
A for the slit width 2a = 150 nm and the screen thickness 
(b = 1000 nm) greater than the wavelength. The figure 
shows the transmission spectra calculated in the study 
and the transmission coefficients calculated by our ana- 
lytical formula and by the numerical evaluation of the 
integral equations 11H14|I for the near- field (z = 150 nm) 
and far-field (z — 1 mm) zones. The transmission spec- 
tra indicates the slit-wave interaction behavior, which 
is similar to those of a Fabry-Perot resonator. The 
transmission resonance peaks, however, have a system- 
atic shift from the Fabry-Perot wavelength Xpp = 2b /n 
(n = 1,2,3...) towards longer wavelengths. We notice 
that the resonance wavelengths determined by our for- 
mulae are in agreement with the model which solves 
the wave equations by using the Rayleigh field expan- 
sion. The peak heights (curves B, C and D), however, 
are different from the results . Our analytical formula 
predicts the great resonant enhancement of a transmitted 



FIG. 5: The transmission as a function of the wavelength 
calculated for the slit width 2a = 150 nm and the screen 
thickness b = 1000 nm. The curve A shows the transmission 
calculated in the study Q]. The curve B shows the trans- 
mission calculated by using the analytical formula 12311 . The 
transmission coefficients calculated by the numerical evalu- 
ation [l3| of the integral equations 111I14H for the far-field 
(z — 1 mm) and near-field (z — 150 nm) zones are shown by 
the curves C and D, respectively. The dotted vertical lines 
correspond the resonant wavelengths given by the analytical 
formula. The respective Fabry-Perot resonances are indicated 
by vertical solid lines. 

wave (curve B), while the Rayleigh expansion model 
indicates attenuation (curve A). Notice that the curves 
B, C and D are practically undistinguishable in the case 
of the narrow slit (2a <C A) . The transmission coefficient 
T calculated by the analutical formula (curve B) and 
the transmission coefficients calculated by the numerical 
evaluation ^| of the integral equations (curves C and 
D) go above unity (i. e., dB) at some resonant wave- 
lengths (see, Fig. 5). Such behaviour of the transmission 
is strange from a point of view of the classical optics. 
Indeed, according to the classical geometrical (A<C2a) or 
Fresnel-Kirchhoff diffraction (A < 2a) optics, the light en- 
ergy impinging on the slit opening only contributes to the 
transmitted energy. The screen completely absorbs the 
light outside the slit aperture. Consequently, the value 
T cannot exceed unity. In the case of the subwavelength 
(A > 2a) slit in a thick (b > A/2) metal screen, the screen 
does not absorb totally the light energy outside the slit 
aperture. At the appropriate resonant conditions, the 
system redistributes the electromagnetic energy in the 
intra-slit region and around the screen, such that more 
energy (T > 1) is effectively transmitted compared to the 
energy impinging on the slit opening. The total energy of 
the system is conserved under the energy redistribution. 
In our plasmon-less model, the redistribution of light en- 
ergy is caused by the boundary conditions imposed on 
the electromagnetic wave at the metal surface. In the 
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plasmon-based models the light energy is redistributed 
due to the boundary conditions for both the electromag- 
netic and electron waves. It can be noted that the trans- 
mission coefficient {T max ~ A/2-7ra at the resonance) in- 
creases with increasing the wavelength. For instance, in 
the experiment |2j the coefficient T ~ 400 was achieved 
in the microwave spectral region. 

Analysis of the denominator of the formulae (|21I23|) 
shows that the transmission will exhibit the maximums 
around the Fabry-Perot wavelengths Xpp — 2b /n (n — 
1,2,3...). The shifts of the resonance wavelengths from 
the values Xpp = 2b /n are caused by the wavelength de- 
pendent term in the denominator of Eq.|531(see, Eas. 121 1 
1221 . The values of the resonant wavelengths and the 
shifts calculated using the analytical formula l|23|) are in 
good agreement with the values calculated using the an- 
alytical formula of the study |f| . The formula ill'.'ill , how- 
ever, shows the great resonant enhancement of a trans- 
mitted wave, while the Rayleigh expansion model 0] pre- 
dicts attenuation. The difference between the Rayleigh 
and Green-function models is attributed to the influence 
of the near-field subwavelength diffraction, which is not 
considered by the models based on the Rayleigh field ex- 
pansion. Indeed, in the Rayleigh-like expansion mod- 
els the fields are found by matching the intra-slit field 
U%{x, z) with the external fields U\(x, z) and U^{x, z) at 
the slit entrance (z — 0) and exit (z = b). The field 
U%{x, z) is presented as a sum of the infinite number of 
waveguide modes, the cavity mode expansion. In the 
most cases, the field U2(x,z) can be well approximated 
by a lowest order waveguide mode. The field U\(x,z) is 
given by a superposition of the infinite number of plane 
waves, the plane wave expansion. The similar plane wave 
expansion is used for the field Us(x,z). In the case of 
A < 2a, the near-field (z — > 0~ or z — > b + ) diffraction 
is weak. Consequently, the fields Ui(x,z) and U^{x,z) 
can be considered as nearly plane waves. Such waves 
are well described by a few terms (m = 0, 1, ... , M) in 
the plane wave expansion. For instance, the classic Airy 
formula for the transmission coefficient is obtained by us- 
ing the zero-order [M = 0) approximation of the fields, 
Ui(x, z) ~ U\{z) and U^{x,z) ~ U^{z). In the case of 
subwavelength (A > 2a) slits, the near-field diffraction is 
strong. The spatial distribution and energy of the fields 
Ui(x,z) and Ua(x,z) are very different from the plane 
waves. To obtain an accurate solution, one should use a 
great number of terms in the plane wave expansion of the 



fields Ui(x, z) and Us(x, z). Unfortunately, the Rayleigh- 
like field expansion models use only a few terms. Thus, 
the fields and transmission coefficient calculated by us- 
ing such models can be considered as a very rough ap- 
proximation. In the case of the Green's function model, 
the rigorous solutions for the fields U\(x, z), ^(x, z) and 
Us{x, z) are found by numerical evaluation of the inte- 
gral equations (|IIII4|) . In the case of the narrow slits, 
we have found an accurate analytical solution of the in- 
tegral equations, which yields the formulae I|I7I23|I for 
the electromagnetic fields and transmission coefficient. 
IV. SUMMARY AND CONCLUSION 

We have presented an analytical model of the reso- 
nantly enhanced transmission of light through a subwave- 
length nm-size slit in a thick metal film. The simple for- 
mulae for the transmitted electromagnetic fields and the 
transmission coefficient were derived by using the narrow- 
slit approximation and the Green's function formalism 
for the solution of Maxwell's equations. The resonance 
wavelengths are in agreement with the semi-analytical 
model 1], which solves the wave equations by using the 
Rayleigh field expansion. Our formulae, however, show 
great resonant enhancement of a transmitted wave, while 
the Rayleigh expansion model predicts attenuation. The 
transparent explanation of the difference between the re- 
sults of the two models was presented. The difference 
is attributed to the near-field subwavelength diffraction, 
which is not considered by the Rayleigh-like expansion 
models. We believe that the presented analytical model 
gains insight into the physics of resonant transmission 
and localization of light by subwavelength nanoapertures 
in metal films. 
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